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Abstract. Observations of binary inspirals with the proposed Laser Interferom- 
eter Space Antenna (LISA) will allow us to place bounds on alternative theories 
of gravity and to study the merger history of massive black holes (MBH). These 
possibilities rely on LISA'S parameter estimation accuracy. We update previous 
studies of parameter estimation for inspiralling compact binaries of MBHs, and 
for inspirals of neutron stars into intermediate-mass black holes, including non- 
precessional spin effects. We work both in Einstein's theory and in alternative 
theories of gravity of the scalar-tensor and massive-graviton types. Inclusion of 
non-precessional spin terms in MBH binaries has little effect on the angular res- 
olution or on distance determination accuracy, but it degrades the estimation of 
intrinsic binary parameters such as chirp mass and reduced mass by between one 
and two orders of magnitude. The bound on the coupling parameter ubd of scalar- 
tensor gravity is significantly reduced by the presence of spin couplings, while the 
reduction in the graviton-mass bound is milder. LISA will measure the luminos- 
ity distance of MBHs to better than ~ 10% out to z ~ 4 for a (10 6 + 10 6 ) M Q 
binary, and out to z ~ 2 for a (10 7 + 10 7 ) Mq binary. The chirp mass of a 
MBH binary can always be determined with excellent accuracy. Ignoring spin 
effects, the reduced mass can be measured within ~ 1 % out to z = 10 and be- 
yond for a (10 6 + 10 6 ) Mq binary, but only out to z ~ 2 for a (10 7 + 10 7 ) M & 
binary. Present-day MBH coalescence rate calculations indicate that most de- 
tectable events should originate at z ~ 2 — 6: at these redshifts LISA can be used 
to measure the two black hole masses and their luminosity distance with suffi- 
cient accuracy to probe the merger history of MBHs. If the low-frequency LISA 
noise can only be trusted down to 10 — 4 Hz, parameter estimation for MBHs (and 
LISA'S ability to perform reliable cosmological observations) will be significantly 
degraded. 
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1. Summary of results 

The Laser Interferometer Space Antenna (LISA) is being designed to detect 
gravitational wave (GW) signals in the frequency band between 10 -4 Hz and 
10 _1 Hz pp. Operating at these low frequencies, LISA can detect, among other 
sources, the inspiral of stellar-mass compact objects - such as neutron stars (NS) 
or black holes (BH) - into intermediate-mass black holes (IMBH) with masses in the 
range 10 2 -10 4 M Q . Another strong gravitational wave source to be observed by LISA 
is the inspiral and merger of massive black holes (MBH) with masses in the range 
10 4 -10 7 M Q . 

Gravitational radiation reaction drives the inspiral of stellar-mass compact objects 
into IMBHs. It also dominates the final stages in the evolution of coalescing MBH 
binaries. The amplitude and phase of the gravitational wave signal carry information 
about binary parameters, such as masses and spins, and about the location and 
distance of the binary. They may also be different in different theories of gravity. 
Therefore LISA can provide important astrophysical information, yield interesting 
tests of fundamental physics, and place bounds on alternative theories of gravity. 
In this paper, we consider, along with standard general relativity, theories of the 
scalar-tensor type (the simplest exemplar being that of Brans and Dicke) and theories 
with an effective mass in the propagation of gravitational waves (which we call 
massive graviton theories, for short). In scalar-tensor theories the phasing evolution 
is modified predominantly by the presence of dipole gravitational radiation reaction 
in the orbital evolution (in general relativity the lowest radiative multipole moment 
is the quadrupole). In massive graviton theories the gravitational wave propagation 
speed depends on wavelength: this generates a distortion in the time of arrival (and in 
the wave phasing) with respect to general relativity, similar to the dispersion of radio 
waves by interstellar plasma. 

In a recent paper |2] (henceforth BBW) we investigated the effect of spin-orbit and 
spin-spin couplings both on the estimation of astrophysical parameters within general 
relativity, and on bounds that can be placed on alternative theories. We restricted 
our analysis to non-precessing spinning binaries, i.e. binaries whose spins are parallel 
(or anti-parallel) to the orbital angular momentum. Our work extended results of 
previous papers [3 EI El El that derived bounds on the graviton mass, on the Brans- 
Dicke parameter wbd and on parameters describing more general scalar-tensor theories 
under the assumption that the compact objects do not carry spin. 

Within Einstein's general relativity, various authors have investigated the 
accuracy with which LISA can determine binary parameters including spin effects 
l^l llUllTTlll2| . These earlier works (except for [T2J) adopted analytical approximations 
to LISA's instrumental noise [5j, augmented by an estimate of white-dwarf confusion 
noise in the low-frequency band. In BBW (and in the present work) we model 
the LISA noise curve by a similar - albeit slightly updated - analytical approximation 
(|14|: see Sec. IIC and Fig. 1 of BBW for details). This noise curve has the advantage 
of being given in analytical form, and reproduces very well the salient features of 
numerical noise curves available online from the LISA Sensitivity Curve Generator |15| , 
a tool sponsored by the LISA International Science Team. 

The central conclusions of our work are as follows. Inclusion of non-precessing 
spin-orbit and spin-spin terms in the gravitational wave phasing generally reduces 
the accuracy with which the parameters of the binary can be estimated. This is 
not surprising, since the parameters are highly correlated, and adding parameters 
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effectively dilutes the available information. Such an effect has already been described 
within Einstein's general relativity in the context of ground-based detectors of the 
LIGO/VIRGO type 117 |. For example, for massive black-hole binaries at 3 Gpc, 
we find that including spin-orbit terms degrades the accuracy in measuring chirp mass 
by factors of order 10, and in measuring the reduced mass parameter by factors of 
order 20 - 100; including spin-spin terms further degrades these accuracies by factors 
of order 3 and 5, respectively. For neutron stars inpiralling into IMBHs with masses 
between 10 3 and 10 4 solar masses, the corresponding reductions are factors of order 20 
and 5 - 30 in chirp mass and reduced mass parameter, respectively, when spin-orbit is 
included, and additional factors of order 4 and 7, respectively, when spin-spin terms 
are included. 

When we consider placing bounds on alternative theories of gravity, for technical 
reasons, we treat only spin-orbit terms. The source of choice to place bounds on the 
coupling parameter wbd of scalar-tensor gravity is the inspiral of a neutron star into 
an IMBH. We first reproduce results of earlier work [B], apart from small differences 
arising from corrected normalization of the LISA noise curve. Then, we include spin- 
orbit effects and find that the bound on ojbd is significantly reduced, by factors of 
order 10 -20. For example for a 1.4 M® neutron star inspiralling into a 400 M Q BH, 
the average bound on wbd from a population of binaries across the sky goes from 
3 x 10 5 to 2 x 10 4 when spin-orbit terms are included. The latter (average) bound 
should be compared with the bound of 4 x 10 4 from Cassini measurements of the 
Shapiro time delay [T%| . 

The effect of including spin on bounding the graviton mass is more modest. In 
this case, the source of choice is the inspiral of MBH binaries. For masses ranging from 
10 5 M Q to 10 7 M Q , the reduction in the bound induced by the inclusion of spin-orbit 
terms is only a factor of 4 to 5. 

We also consider the effect of spin terms on the angular and distance resolution 
of LISA. We find that spin couplings have a mild effect on the angular resolution, on 
the distance and, as a consequence, on the redshift determination for MBH binaries. 
By contrast, for stellar mass objects inspiralling into IMBHs, neither distance nor 
location on the sky is very well determined. Furthermore, for NS-IMBH binaries the 
angular resolution is somewhat dependent on the inclusion of spin terms. Here we 
show that the different spin dependence of the angular resolution in the NS-IMBH 
and MBH-MBH cases can be traced back to a different behavior of the correlations 
between the chirp mass (leading the gravitational wave evolution of the binary) and 
the angular variables describing LISA's position in the sky. 

LISA's low-frequency sensitivity affects the accuracy of estimating parameters 
for MBH binaries, as we showed in BBW (similar investigations can be found in 19; 
and 20 ). As default low-frequency cutoff we choose 10 -5 Hz. However, below 
10 -4 Hz, the noise characteristics of LISA are uncertain, and if we set LISA's low- 
frequency cutoff at 10~ 4 Hz, the accuracy of estimating both extrinsic parameters 
such as distance and sky position and intrinsic parameters such as chirp mass and 
reduced mass, as well as the graviton mass, can be significantly lower, especially for 
higher-mass systems. 

LISA can observe MBH binaries with large SNR out to large values of the 
cosmological redshift. If the corresponding mass and distance determinations are 
accurate enough, LISA will be a useful tool to study BH formation in the early 
Universe. Using Monte Carlo simulations we find that LISA can provide accurate 
distance determinations out to redshift z ~ 2 for source masses ~ 10 7 M Q , and out 
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to z ~ 4 for source masses ~ 10 6 Af Q . Mass determinations strongly depend on an 
accurate treatment of spin effects. The chirp mass of a MBH binary can always be 
determined with excellent accuracy. Ignoring spin effects, the reduced mass can be 
measured within ~ 1 % out to z = 10 and beyond for a (10 6 + 10 6 ) M© binary, 
but only out to z ~ 2 for a (10 7 + 10 7 ) M binary. Present-day MBH coalescence 
rate calculations indicate that most detectable events should originate at z ~ 2 — 6 
[211 1221 I23| : at these redshifts LISA can be used to measure the two BH masses 
and their luminosity distance with sufficient accuracy to probe the merger history of 
MBHs. 

The paper is organized as follows. In Sec. [21 we briefly summarize our procedure 
to study parameter estimation with LISA. All calculations presented here adopt the 
"more realistic" of the two models described in BBW: we take into account the orbital 
motion of the detector and use Monte Carlo simulations of a population of sources 
across the sky. In Sec. [21 and 0] we extend and clarify some results presented in BBW. 
Sec. [31shows the dependence of parameter estimation and tests of alternative theories 
of gravity on the binary mass (BBW only considered a few representative values of the 
masses). Sec. 0] discusses the cosmological reach of LISA and its ability to probe the 
merger history of MBHs, extending and complementing the discussion in Sec. IIIC of 
BBW. 

2. Gravitational wave detection and parameter estimation 

We assume, as in that two independent Michelson outputs can be constructed 
from the readouts of the three LISA arms if the noise in the arms is totally symmetric 
(but see [21] for discussion of the most general combinations of LISA outputs, and 
their sensitivities). Since the modulation of the measured signal due to LISA's motion 
occurs on timescales ~ 1 year (much larger than the binary's orbital period), the 
Fourier transform of the signal can be evaluated in the stationary phase approximation, 
with the result 

gi('W)-¥'i»,a 1 (t)-¥>D(t)) , (1) 

Here / is the frequency of the gravitational waves; M. = r) 3 / 5 M is the "chirp" mass, 
with M = mi + TO2 and 77 = m\milM 2 \ is the luminosity distance to the source; 
A a (t) = {[1 + (L ■ n) 2 ] 2 F+ 2 + 4(L • n) 2 F Q X 2 } 1/2 , where L is the orbital angular 
momentum unit vector, n is a unit vector in the direction of the source on the sky, 
and the quantities F^ ,x are the pattern functions for the two equivalent Michelson 
detectors. The time variable t = t(f) is given by an expansion to 2PN order, including 
also the Brans-Dicke parameter and the graviton-mass term (see BBW and [S] for 
details); <p p , a (t) is the waveform polarization phase and <fD(t) the Doppler phase. 

We have adopted the standard "restricted post-Newtonian approximation" for 
the waveform, in which the amplitude is expressed to the leading order in a post- 
Newtonian expansion (an expansion for slow-motion, weak-field systems in powers of 
v ~ (M/r) 1 ' 2 ~ (7T.M/) 1 / 3 ), while the phasing ip(f), to which laser interferometers 
are most sensitive, is expressed to the highest post-Newtonian (PN) order reasonable 
for the problem at hand. For binaries with spins aligned (or anti-aligned) and normal 
to the orbital plane, this is a valid approximation because the amplitude varies slowly 
(on a radiation reaction timescale) compared to the orbital period. But when the 
spins are not aligned, modulations of the amplitude on a precession timescale must be 
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included. Preliminary studies show that precession effectively decorrelates parameters, 
improving parameter estimation in a significant way The study of precessional 
modulations is beyond the scope of this work. 

The phasing function ip(f) is known for point masses up to 3.5 PN order |25II26| . 
But spin terms are known only up to 2PN order, so to be reasonably consistent, we 
include in the phasing point-mass terms only up to this same 2PN order. The needed 
expression for the phasing is 

1>(f) = 2tt ft c ~ (f> c (2) 

^2 



128 tt 2 DM , .,,,,2/3 /3715 55 \ _ 2/5 , ..,, 2/3 

- - AfTTT^y (7rM/) 1 + [-756 + yV 11 1 ^ Mf) / 

- Mirrj-V 5 (irMf) + 4/3 rT 3/5 {^Mf) - Vdan^ 15 (irMf) 4/3 
/ 15293365 27145 3085 2 \ ... . ,^, 4/3 ] 

where t c and <p c are the time and phase at coalescence. The second term inside the 
braces is the contribution of dipole gravitational radiation in Brans-Dicke theory: it 
depends on the difference of the scalar charges of the two bodies af ~ af/2uiBD 
(* = 1,2) §. 

For black holes, a, = 0, therefore we need "mixed" binaries (eg. NS-IMBH 
binaries) to place bounds on scalar-tensor theories. The third term in the braces is 
the effect of a massive graviton, which alters the arrival time of waves of a given 
frequency, depending on the size of the graviton Compton wavelength X g and on a 
distance quantity D (see BBW or 0] for the relation between D and the luminosity 
distance Dl)- The remaining terms in the braces are the standard general relativistic, 
post-Newtonian terms, including spin effects. 

The quantities (3 and a represent spin-orbit and spin-spin contributions to the 
phasing, given by 



1 2 
^=12 ^> 

i=l 



113^+75, 



L • Si , (3) 



<r = ^XiX2 (-247Si • S 2 + 721L • S X L • S 2 ) , (4) 

where Si and L are unit vectors in the direction of the spins and of the orbital 
angular momentum, respectively, and Si = XifTifSj. For BHs, the dimensionless spin 
parameters Xi must be smaller than unity, while for neutron stars, they are generally 
much smaller than unity. It follows that \/3\ < 9.4 and |er| < 2.5. 

We use the by now standard machinery of parameter estimation in matched 
filtering for gravitational wave detection that has been developed by a number of 
authors ^3 1281 1291 13D] . Given a set of parameters 9 a characterizing the waveform, 
one defines the "Fisher matrix" r ao with components given by 



r„ 



/ dh 


dh_\ 
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'dh* dh 




de b J 


~ 2 Jo S n (f) 





dh* dh 
~d¥~dlF 



df. (5) 



§ Note that the quantity Oi = 1 — 2sj, where s; is the so-called sensitivity of the body 0. In 
calculations involving NSs we adopt the value ajjg = 0.6, a typical value obtained from calculations 
of NS structure with realistic equations of state in scalar- tensor theories 1271 . 
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An estimate of the rms error, A9 a , in measuring the parameter 9 a can then be 
calculated, in the limit of large SNR, by taking the square root of the diagonal elements 
of the inverse of the Fisher matrix, 

A9 a = V^, £ = r" 1 . (6) 

The correlation coefficients between two parameters 9 a and 9 b are given by 

c ab = S a VVS aa E &f) . (7) 

The assumption of large SNR is well justified for MBH binaries, which can be observed 
with high SNR out to very large distances. 

To determine how accurately LISA can measure source locations and luminosity 
distances we perform Monte Carlo simulations using a population of sources across the 
sky. We consider in detail two classes of systems. The first are NS-IMBH binaries with 
a NS mass Mns = 1-4 M© and a BH mass in the range [400— 10 4 ] M©, a typical target 
system used to place bounds on the BD parameter. According to current estimates, 
detection rates for 10 - 10 2 M© BHs spiralling into IMBHs of mass 10 2 - 10 3 M 
are expected to be very low |31j . In order to apply the Fisher matrix formalism in a 
consistent way, we assume that these NS-IMBH binaries, rare though they may be, 
are detected with a single-detector SNR pj = 10. The second arc MBH binaries. In 
this case we first consider equal-mass binaries with total measured mass in the range 
[2 x 10 4 — 2 x 10 7 ] Mq at distance Dl = 3 Gpc - typical systems to put bounds 
on massive graviton theories (Sec.|2J). Then we study parameter estimation for MBH 
systems located at different cosmological redshifts and having mass (10 6 + 10 6 ) M© 
and (10 7 + 10 7 ) M© as measured in the source rest frame (Sec.QJ. 

For each of these systems we distribute 10 4 sources over sky position and 
orientation, described, in the solar system barycentric frame, by angles (9s,4>s) and 
(9~l,4>l), respectively. We randomly generate the angles (f>s, 4>l in the range [0, 27r] 
and ps = cos 9s, — cos9~l in the range [—1,1] (see BBW for details). When the 
waveform contains a large number of highly correlated parameters, computing the 
inverse of the Fisher matrix can be numerically difficult. The method we used to 
check the robustness of our results is described in Appendix B of BBW. 

3. Parameter estimation and bounds on alternative theories of gravity 

In BBW we carried out a set of Monte Carlo simulations for a MBH binary of mass 
(10 6 + 10 6 ) Mq at fiducial distance Dl — 3 Gpc. We showed that our results are 
compatible with results obtained using a simpler approach in which we don't take 
into account LISA'S orbital motion and we average over pattern functions. 

For this paper, we performed a more extensive set of Monte Carlo simulations, 
spanning a wider range of binary masses. Results from these simulations are shown 
in Fig. n] We consider nonspinning equal-mass MBHs at fiducial distance Dl = 
3 Gpc, and we investigate the dependence of mass, distance and angular resolution 
determinations on the binary's total mass. From the plots on the first row we see 
that distance and angular resolution determinations are largely independent of the 
inclusion of spin terms (and even of the massive-graviton terms) in the phase. Mass 
determinations are more sensitive to the inclusion of spin terms and to the theory of 
gravity we arc considering. Distance errors and angular resolution generally decrease 
with the binary's mass. On the contrary, mass errors have a relative minimum in the 
given mass range. This is because the measurement of mass parameters depends on 
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Figure 1. Top: average errors on Cg (dashed lines) and Dj, (solid lines). Bottom: 
average errors on M (dashed lines) and /i (solid lines). Left panel refers to GR, 
right panel to massive graviton theories; errors are given as a function of the total 
mass for equal mass MBH binaries. We assume the LISA noise curve can be 
trusted down to /j ow = Hz. Black lines are computed omitting spin terms, 

red lines include spin-orbit terms, blue lines include both spin-orbit and spin-spin 
terms. 



an accurate modeling of the phase of the waveform: when the MBH mass is large the 
binary's frequency becomes low, the binary spends a significant fraction of time out 
of the LISA sensitivity band, and the accuracy of mass determinations decreases. 

We also explored the dependence of parameter estimation on the signal's ending 
frequency. Instead of using the Schwarzschild Innermost Stable Circular Orbit (ISCO) 
as a conventional ending frequency of the signal, we looked at the effect of truncating 
the signal at the Minimum Energy Circular Orbit (MECO; see Sec. IIB of [SJj for 
a discussion), consistently computed at 2PN and assuming dynamically small spins 
(/3 = a = 0). We found that using the MECO instead of the ISCO has a completely 
negligible effect on distance and angular position determinations. The effect on mass 
determinations is completely negligible for total masses M smaller than ~ 10 6 M Q ; 
even for values of M larger than this, using the MECO instead of the ISCO only 
results in slightly smaller (by factors ~ 2 or so) errors. However, the effect of using 
the MECO could be non-negligible if the spins are large. 

In Fig. |21 we look at bounds to be placed on i) the graviton Compton wavelength 
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Figure 2. Left: average bounds on the graviton wavelength as a function of 
total mass for equal mass MBH binaries. Right: bound on the BD parameter 
from a NS-IMBH binary as a function of the IMBH mass. The horizontal blue 
line corresponds to the Cassini bound 1181 . Black lines are computed omitting 
spin terms; red lines include spin-orbit terms as well. The dashed line shows that 
the bound on A 9 is reduced if the LISA noise curve can only be trusted down to 
/low = 10- 4 Hz. 



X g , using a sample of equal-mass MBHs at fiducial distance Dl = 3 Gpc (left panel), 
and ii) the Brans-Dicke parameter u>bd, using a sample of NS-IMBH binaries with 
single-detector signal-to-noise ratio equal to 10, as a function of the IMBH mass 
(right panel) . These results confirm qualitatively (and improve quantitatively) results 
previously obtained by Will and collaborators. Those authors adopted a somewhat 
cruder LISA model, that does not take into account LISA's orbital motion O0J|SJE1- 
Fig. |2 shows that understanding the low-frequency LISA noise is important to set 
bounds on the graviton mass using high-mass MBH binaries. We will see in Sec. 0] 
that an understanding of the low frequency noise is also important to use LISA in a 
cosmological context. 

In BBW we found that LISA's angular resolution degrades when we include spin 
terms in the NS-IMBH case, but it is essentially unaffected in the MBH case. To 
explore this further, in Fig. [3] we show a scatter plot (out of 10 4 binaries) of the 
correlation c^ sM between the azimuthal angle 4>s (describing the source position in the 
sky with respect to the Solar System barycenter) and the chirp mass M. (dominating 
the phase evolution of the gravitational waveform), as a function of 4>s, in the two 
cases. The correlation has a "clean" periodic behavior in the NS-IMBH case; in the 
MBH-MBH case Cj, sj ^ is not a clean periodic function of the binary's azimuthal angle 
4>s, showing some beating phenomena. An important feature is that, when we include 
spin terms, the maximum correlation between mass and angular position clearly grows 
in the NS-IMBH case, but not in the MBH-MBH case. This is consistent with the 
fact that the angular resolution depends on spin couplings in the NS-IMBH case, but 
it doesn't in the MBH-MBH case. 
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Figure 3. Scatter plot of the correlations c<f, s jvi{4>s) m GR. From left to right 
we observe how correlations change when we add spins. The top row corresponds 
to a (10 6 + 10 6 ) Mq MBH binary at Dj, = 3 Gpc; the bottom row refers to a 
(1.4 + 10 3 ) Mq NS-IMBH system observed with single-detector SNR equal to 10. 



4. Probing the merger history of supermassive black holes with LISA 

Astrophysical estimates of MBH coalescence rates depend on many complex physical 
phenomena. The study of MBH binary evolution in galaxy mergers is a very active 
field of research [jJ5J • A class of models predicts that the LISA mission should be able 
to detect at least ~ 10 MBH coalescence events out to large values of the redshift, 
and that most of these events would originate at z ~ 2 — 6 E5J . Research 

in this field is very active and the estimates are still quite uncertain: a partial list of 
references includes |51 1551 1551 1771 1551 155] . 

To use LISA as a cosmological probe it would be desirable to have an accurate 
measurement of the binary parameters, and not just a simple detection. In this 
context the redshift z can be large, and the distinction between observed masses 
and masses as measured in the source rest frame - given by the simple rescaling 
-^observed = (1 + z)M sourcc , where z is the cosmological redshift - is important||. LISA 
can only measure redshifted combinations of the intrinsic source parameters (masses 



Up to this point we denoted by M and M the observed chirp and total masses. 
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and spins), so it cannot measure the redshift z to the source. But if the MBH masses, 
luminosity distances and angular resolution are determined with sufficient accuracy, 
LISA can still provide information on the growth of structures at high redshift. If 
cosmological parameters are known, the relation between luminosity distance and 
redshift Dl(z, Hq, SIm, ^a) (where Hq = 72 km s _1 Mpc" 1 is Hubble's constant and 
we assume a zero-spatial-curvature Universe with Q K = 0, 57a + &m = 1 5 according 
to the present observational estimates) can be inverted to yield z(Dl, Hq, Qm, ^a)- 
Then LISA measurements of the luminosity distance can be used to obtain intrinsic 
BH masses as a function of redshift, thus constraining hierarchical merger scenarios UJ. 
Alternatively: if we can obtain the binary's redshift by some other means, e.g. from an 
electromagnetic counterpart, then LISA measurements of Dl can be used to improve 
our knowledge of the cosmological parameters |4()l 1411 IB?] . 

These exciting applications depend, of course, on LISA's measurement accuracy 
at large redshifts. Here, as in BBW, we look at the redshift dependence of 
measurement errors for two representative MBH binaries having masses (10 6 + 10 6 )Af Q 
and (1O 7 + 1O 7 )M0 as measured in the source rest frame. We also assume that the LISA 
noise can be extrapolated down to /i ow = 10~ 5 Hz; more conservative assumptions on 
/low could significantly affect our conclusions (see |19l 1201 12*|). We compute the errors 
performing Monte Carlo simulations of 10 4 binaries for different values of the redshift 
and then averaging over all binaries. 

As discussed in Sec. OH distance determination and angular resolution for MBHs 
are essentially independent of the inclusion of spin terms. The (average) relative error 
on D L for the (10 6 + 10 6 )M Q system is ~ 2% at z = 1, ~ 5% at z = 2 and - 11% 
at z = 4. This reduction in accuracy is due to the fact that the signal spends less 
and less time in band as the redshift is increased. We only consider values of the 
redshift such that the binary spends at least one month in band before coalescing: 
for the case (10 7 + 10 7 )M©, this corresponds to z ~ 4. The distance determination 
error for this high-mass binary grows quite rapidly, being ~ 2% at z = 1, ~ 6% at 
z = 2 and ~ 21% at z = 4. LISA's angular resolution is rather poor even at small 
redshifts, and it rapidly degrades for sources located farther away, the degradation 
being more pronounced for higher-mass binaries. Better distance determinations can 
be obtained if we are lucky enough to locate the source in the sky by some other 
means: for example, associating the gravitational wave event with an electromagnetic 
counterpart. In this case angles and distance would be decorrelated, allowing order of 
magnitude improvements in the determination of Dl ^Hj- For (10 6 + 10 6 )M Q general 
relativistic nonspinning binaries, a least-square fit of mass and distance errors in the 
interval z € [1, 10] yields: 

AM /M = (-0.1148 + 0.7236z + 0.5738z 2 ) x 10~ 5 , 

Afi/fj, = (-0.61431 + 1.9018^ + 0.43721z 2 ) x 10~ 4 , 

AD L /D L = (-0.65651 + 2.6935z + 0.061595z 2 ) x 10~ 2 . 

It is important to stress that in our discussion above we refer to average errors from a 
population of 10 4 sources randomly distributed across the sky. In a single "lucky" 
detection, parameter estimation errors could be much smaller than the average. 
The most useful quantity to measure LISA's cosmological capabilities is perhaps the 
percentage of binaries whose parameters can be measured within some given accuracy, 
as a function of z. This percentage is shown in Fig.0J 

Upper limits on measurement errors of the various parameters have been chosen 
in order to answer the following questions. Out to which value of z can we determine 
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Figure 4. Percentage of binaries for which: Sfi = A/x/yU < 1% (black circles), 
SM = AM/M < 0.1% (red squares; 8M = AM/M < 1% in all cases 
considered), 8D L = AD L /D L < 10% (blue diamonds), SD L = AD L /D L < 5% 
(blue diamonds, dashed lines), AQ$ < 10~ 4 (green triangles). On the left we 
consider a BH-BH binary of mass (10 6 + 10 6 ) Mq, on the right a BH-BH binary 
of mass (10 7 + 10 7 ) Mq. We assume / low = 10 -5 Hz. 



the chirp mass to better than 0.1 %, and the reduced mass to better than 1 %? 
Out to which z does LISA beat the present uncertainty (~ 10%) on cosmological 
parameters, allowing us to determine the luminosity distance to the source? Which 
fraction of binaries can be observed with an angular resolution error AQ$ < 10 -4 , 
roughly corresponding to the angular diameter of the Moon (and the Sun) as seen 
from the Earth? 

The answer to these questions depends of course on the binary's mass and on 
the low- frequency sensitivity of LISA. Following BBW, in our calculations we assume 
that the noise curve can be extrapolated down to f\ ow = 10~ 5 Hz and we consider two 
representative binaries with mass (10 6 + 10 6 ) Mq and (10 7 + 10 7 ) Mq, respectively. 

The main conclusions to be drawn from Fig.0]are the following, (i) The fractional 
error on the chirp mass can be better than 0.1 % out to z = 10 and beyond for a 
(10 6 + 10 6 ) Mq binary, but only out to z ~ 2 for a (10 7 + 10 7 ) M© binary, (ii) 
In the absence of spins, the fractional error on the reduced mass can be better than 
~ 1 % out to z = 10 and beyond for a (10 6 + 10 6 ) Mq binary, but only out to z ~ 2 
for a (10 7 + 10 7 ) Mq binary, (iii) For ~ 80 % of observed binaries, LISA's distance 
determination will be better than ~ 10% (beating present errors on cosmological 
parameters) out to z ~ 4 for a (10 6 + 10 6 ) Mq binary, and out to z ~ 2 for a 
(10 7 + 10 7 ) Mq binary, (iv) LISA's angular resolution is comparatively quite poor: 
only a small fraction of binaries can be observed with Af^s < 10~ 4 at any redshift. 

Unimportant as they are for distance determination and angular resolution, spin 
effects have a dramatic impact on mass measurement accuracy. In BBW we showed 
that when we include spin-orbit and spin-spin terms, the error on the chirp mass is still 
pretty small: it only becomes larger than a few percent for massive (M ~2x 10 7 M Q ) 
binaries located at z > 2. However, non-precessional spin effects seriously limit our 
ability to measure the reduced mass. When we include both spin-orbit and spin- 
spin terms, the reduced mass error becomes larger than 5 % at all redshifts even for 
"low mass", (10 6 -I- 10 6 ) M Q binaries. Our conclusions on mass determinations for 
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spinning binaries could change when we consistently include ^recessional spin effects. 
These effects can induce modulations in the waveform, possibly improving the mass 
measurements in a significant way |llj . Therefore the study of precession is very 
important to assess LISA'S ability to measure MBH masses in galactic mergers; we 
plan to return to this study in the future. 
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